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Abstract We investigate the radial color gradients of galactic disks using 
a sample of ~ 20, 000 face-on spiral galaxies selected from the fourth data 
release of the Sloan Digital Sky Survey (SDSS-DR4). We combine galax- 
ies with similar concentration, size and luminosity to construct composite 
galaxies, and then measure their color profiles by stacking the azimuthally 
averaged radial color profiles of all the member galaxies. Except for the 
smallest galaxies (i?5o < 3 kpc), almost all galaxies show negative disk 
color gradients with mean g — r gradient Ggr = —0.006 mag kpc^^ and 
r — z gradient Grz = —0.018 mag kpc^^. The disk color gradients are in- 
dependent of the morphological types of galaxies and strongly dependent 
on the disk surface brightness /i^, with lower surface brightness galactic 
disks having steeper color gradients. We quantify the intrinsic correlation 
between color gradients and surface brightness as Gg,- = —0.011/2^+0.233 
and Grz = —0.015/1^ + 0.324. These quantified correlations provide tight 
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observational constraints on the formation and evolution models of spiral 
galaxies. 



Key words: galaxies: spiral - galaxies: statistics - galaxies: evolution - 
galaxies: fundamental parameters (color gradient, luminosity, radius, sur- 
face brightness) 



1 INTRODUCTION 



Radial color gradient of a spiral galaxy freezes its fossil information of the star for- 
mation history along the disk and therefore provides a strong constraint on the mod- 
els of disk formation and evolution. Earlier studies revealed that most of the spiral 
galaxies show negative color gradi e nts, i.e. the color is gradually bluer outwards (e.g. 



Kim & Ann 



V negative color gradi e nts, i.e. tne color is gradually bluer outwards (e.g. 
199({ Ide Jonglll996bl: Isell & de Jonsl llOod: loadotti & dos AnjoslboOl 



Moth & Elston 



Kim & Ann 



200: 



i 



MacArthur et al. 



20041 : 



Taylor et al 



20051). 



(119901) studied the radial color distributions of a sample of 103 face- 
on galaxies with different morphological types. They concluded that about 40 percent 
of elliptical and lenticular galaxies and about 75 percent of intermediate type spiral 
galaxies show clear negative color gradients, while the latest type spiral and irregular 
galaxies have positive gradients or no gradient at all. For spiral galaxies, they found a 
weak correlation between the color gradient and absolute magnitude, with more lumi- 
nous spirals having steeper gradients. However, in their work, the color gradient was 
measured along the whole image of each galaxy, where the bulge and disk components 
were not decomposed. 

It is well established that a spi ral galaxy can be decomposed into a concen 



trated bulge plus an exp onential disk (|de Vaucouleurs 



1959 



Simien & MichardL 



1990 



Capaccioli et al.Lll992h . Therefore, the negative color gradient of a spiral galaxy might 
be partially originated from the configuration of these two components. To remove this 
effect and measure the color gr adients of these two components separately, a bulge- 
disk decomposition is required, de Jong (|l996a b ^ made the bulge-disk decomposi- 
tion for a sample of 86 face-on spiral galaxies and investigated radial color distribu- 
tions of disk components in optical and near-infrared bands. He found that the color 
of disk component also becomes bluer outwards, but no significant correlations were 
found between the color gradients and other structural parameters of galaxies, such 
as the color, inclination, morphological type, central surface brightness, scale-length, 
integrated magnitude and so on. However, the decomposition is dependent on the mod- 
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els use d. Different models used could give different results. Alternatively. iTaylor et al 



(l2005h investigated the color gradients of the bulge dominated inner parts and the disk 
dominated outer parts of 142 galaxies by simply separating each galaxy at the effec- 
tive radius Re. For the disk parts, they found that the early and intermediate type spiral 
galaxies show weak negative color gradients, while the late type spiral and irregular 
galaxies have positive color gradients. 

Although the color gradient has been investigated extensively, its quantitative def- 
initions are still diverse in different studies. There are mainly two branches of such 
definitions. One is the color variation AC I as a functio n of scaled radiu s (e.g., scale- 
length -Rrf, half-light radius -R.^n), either in li near space (ISegalovia. 



1975 



1988 



2003 



Kim & Ann, 



IWu et al 



1990; 



2005 



Vader et al.. 



Taylor et al.1. |2005|) or in logarithmic space (ITamura & Ohta . 



iLi et al, 



20051). This kind of definition implies that galaxies 



with different sizes have similarities in the structure of color distribution. The other 
bran ch is the AC/ measured along the physical radius, e.g. the color variation pe r 



kpc (iPrugniel & Heraudeaul 1 19981: iGadotti & dos Anjos 



2001 



Moth & Elston 



2002h . 



which could also be presented in either linear or logarithmic space. Additionally, if 
one assumes that the galactic disk in different wavebands could all be well modeled 
by exponential profiles, the change of scale-len gth as a function of wavebands could 



also be used t o parameterize the color gradient (|Peletier et al.L 1 19941 : Ide JongL 1 1996b 



CunowL l2004h . Despite the diversity of the definitions, the color gradient is a essential 



parameter that quantify the systematical change of color profile of galaxy. 

To derive the color profile of a galaxy, detailed photometric analyses in different 
wavebands are required. Moreover, to get robust conclusions on the correlations be- 
tween color gradients and other structural parameters, a large and homogenous sam- 
ple of galaxies is needed. Up to now, the number of galaxies with detailed measure- 
ment of color profile is qui t e limited (IGadotti & dos AniosL 1200 ll : IPohlen. & Trujillo , 



2006: Munoz-Mareos et al, 



2007 



AzzoUini et al 



York et al. 



2000l) has surveyed about a quarter of the whole 



2008b . On the other hand, the Sloan 



Digital Sky Survey (SDSS, 
sky in 5 broad wave bands (u, q, r, i, z) and provided a large catalogue of photomet- 
ric data of galaxies (jAdelman-McCarthy et al.L 120071) . In the released catalogue, the 
surface brightness profile of each galaxy h as been measured in a series of fixed con- 
centric annuli, known as prof Meai^ (see Stoughton et al. 2002 for detail, hereafter 
EDR). Due to the short exposure time (about 54 second) of SDSS image, the signal- 
to-noise ratio (S/N) of prof Mean data is limited and it is suggested to be only used in 



' Before calculating the surface brightness profile prof Mean, the fore gro und and background ob - 
jects have been deblendedby the FRAMES pipeline, see lLupton et aLloOOll) and lStoughton et alJ(l2002h . 
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the cumulative profile for individual galaxies in EDR papeo However, the S/N could 
be greatly improved by stacking the profiles (prof Mean) of a large number of simi- 
lar galaxies. The large amount of galaxies in SDSS provides an opportunity to make 
subsamples with enough member galaxies for such stacking. Comparing to the direct 
image stacking method with contamin ated sources carefully masked (e.g. 



Zibetti et al. 



20041 : Ide Jongll2008l : iHathi et al.ll2008|) . this kind of catalog data stacking is more vul- 



nerable to background co ntaminations, but is much simpler and mo re suitable for the 



analysis of large data sets (jAzzoUini et al 



2008 



Bakos et al. 



2008). 



In this study, we use a large sample of face-on spiral galaxies (~ 20,000) selected 
from SDSS-DR4 (|Adelman-McCarthy et al.L 120061) to investigate their disk color gra- 
dients statistically. We divide these spiral galaxies into subsamples with similar con- 
centrations, sizes and luminosities, and then construct composite galaxies by combin- 
ing the color profiles of the member galaxies in each subsample. After measuring the 
color gradients of disk components of these composite galaxies, we further explore 
possible correlations between color gradients and structural parameters. The main aim 
of this work is to find out which parameter is primarily correlated with the color gra- 
dient and discuss its physical implications. 

The structure of this paper is as follows. In Section[2l we describe our sample selec- 
tion and the building of composite galaxies. The color gradients and their correlations 
with other structural parameters are presented in Section [3l Finally, we give a brief 
summary and make discussions in Section|4l 



2 DATA 



2.1 Sampling 



The galaxies used in this paper are drawn from the main galaxy sample of the SDSS- 
DR4, which is a spectroscopic sample complete to r-band Petrosian magnitudes rp < 
17.77. We use galaxies with rp < 17.5 and redshift z > 0.01. The magnitude limit 
excludes the faint objects, ensuring the st ructure of the galaxies being nicely resolved 
by SDSS images (|Desroches et al.l . l2007b . The low redshift limit eliminates the effect 
of the peculiar velocity on the distance estimation. 



^ iBlanton et al.l ( l2003bl) has used prof Mean and successfully fitted the surface brightness profile 
with a Sersic model for each individual galaxy. 
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Fig. 1 The histograms of fracDeV [left panel (a)] and axis ratio b/a [right 
panel (b)]. The dotted histograms show the distributions of the initial sample 
with rp < 17.5 and z > 0.01 while the shadowed regions represent the distri- 
butions of final sample of 25,517 face-on spiral galaxies (rp < 17.5, z > 0.01, 
fracDeV < 0.5 and b/a > 0.75). The solid histogram in panel (b) shows 
the b/a distribution of selected spiral galaxies (fracDeV < 0.5), see text for 
detail. 



To select spiral galaxies, we adopt the criterion fracDeV < 0.5 (in r-ban(£|). The 
parameter fracDeV indicates the likelihood of the surface brightness profile that can 
be best modeled by a de Vaucouleurs profile ( Abazaiian et all 



wide 



2006 



V used as a galaxy r norphology indicator (e 



ShaoetaL 



20071 : 



Padilla & Strauss 



2008 



2004 ) and has been 



Bemardi et 



al 



Zhong et alj 



2005 



20081 



i 



Chang et al 



Strateva et al 



(|200l|) have shown that the criterion fracDeV < 0.5 can separate late type galax- 
ies from early type ones with a reliability ~ 90 percent. We show the histogram of 
fracDeV of sample galaxies in the left panel of Fig. 1 . As we can see, for the galaxies 
with fracDeV < 0.5, most of them peak at fracDeV ~ 0. 

Aiming at the study of color gradients, we further restrict our spirals to be only face- 
on galaxies. This restriction is based on two considerations. One is that the internal dust 
attenuation in the face-on case is minimized. The other is that the prof Mean, which 
is measured in circular apertures, is more reliable for face-on galaxies. We adopt the 
criterion b/a > 0.75 to make such face-on selection, where 6/a is the axis ratio derived 



^ In this paper, all the photometric parameters are in r-band as default unless we specify the other 
band with a subscript. 
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from the exponential profile fitting and given as ab_exp in SDSS catalogj (see EDR 
paper). We show the histograms of h/a of our sample galaxies in the right panel of 
Fig.l. The dotted and solid histograms represent the distributions of galaxies before 
and after the morphological selection (fracDeV < 0.5), respectively. As expected, 
the b /a shows a flat distribution and paramet erizes the inclination of spira l disk very 



well (Lambas et al.L 1 19921 : IShao et al.l . l2007h . According to the study of IShao et al 



(|2007h . the axis ratio 6/a ~ 0.75 corresponds to inclination angle i ~ 41° (see fig. 3 
of their paper). 

Finally, we select a magnitude limit sample of 25,517 face-on spiral galaxies, 
whose distributions of fracDeV and b/a are shown as shadow regions in Fig.[IJ 

2.2 Binning and composing galaxies 

In this section, we describe how to bin sample galaxies and construct composite galax- 
ies. To construct a composite galaxy, the member galaxies must have similar physical 
properties. We assume that, if a set of galaxies have similar concentrations (morpho- 
logical type), physical sizes (angular momenturr@) and luminosities (stellar mass), their 
other physical properties are also similar, including the radial color profiles. 

We divide the 25,517 sample galaxies into subsamples in small volumes of concen- 
tration c, physical half-light radius R^q (in unit of kpc) and absolute Petrosian magni- 
tude M. The bin widths of c, R^q and M are 0.1, 1.0 kpc and 0.5 mag, respectively. The 
concentration is defined as c = Rgo/R^o, where Rgo and -R50 are aperture radii con- 
taining 90 and 50 percent of Petrosian flux (see EDR paper), respectively. The physical 
size and absolute magnitude are calculated under the cosmology contex with = 0.3, 
l^A = 0.7 and Hq = 70 km s^^Mpc-^ 

After dividing of galaxies into subsamples, we measure the dispersion of the colors 
of all the member galaxies in each subset. This dispersion is very small, with typical 
value ~ 0.1 mag (comparable to the bin width of the absolute magnitude), which 
confirms our assumption that all the physical properties of member galaxies are similar 
when c, R^m and M are constrained. 

Since the member galaxies of each subset are assumed to have similar color profiles 
but located at different redshifts, their color profiles (measured from prof Mean in a 
series of fixed angular sizes) could be treated as the color measurement of 'one' galaxy 



^ Both fracDeV and ab_exp are point spread function corrected fitting parameters in the SDSS 
catalog. 
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At given magnitude, the galaxies with larger angular momentum (parameterized by larger spin 



parameter A) tend to have larger sizes (IMo. Mao & Whitelll998tlSver. Mao & Moll 19991) 
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at different physical radii. Therefore, the color profile of this composite galaxy could 
be constructed by combining color pro files of all member galaxies in a subset (e.g., 



AzzoUini et al. 



2008 



Bakos et al. 



2008). 



Before make the composite color profiles, we have made a global K-correction 
on the color p r ofile o f each member galaxy usin g the routine kcorrect-v.4. 1 .4 



(|Blanton et al.L I2003al: iBlanton & RoweisL 120071) . In principle, for each member 



galaxy, if color gradient exists, the K-corrections at different radii will be different. 
This second-order effect has been tested by using the prof Mean data of 5 SDSS 
wavebands to make K-corrections for each photometric annulus independently. We 
found the difference between the radius dependent K-correction and the global one is 
less than 0.001 mag. 

Fig. [2] shows g — r and r — z color profiles of an example composite galaxy which 
comprises 166 member galaxies within the volume of 1.9 < c < 2.0, 7.0 < R^q < 
8.0 kpc and —21.5 < M < —21.0 mag. In this figure, each small dot denotes the color 
in a prof Mean annulus of one member galaxy. These small dots are then divided into 
30 bins, with 15 bins inside R^o and the others outside. Each inside bin (R < R^q) 
contains the same number of data points (average number of all small dots inside -R50), 
so do the outside bins (R > R50). We average the colors of small dots in each bin and 
print them with large circles, with open and filled ones representing inside and outside 
-R50 respectively. As we can see, by averaging a group of member galaxies, the color 
profile of the composite galaxy tends to be smooth in the sense that the S/N of color 
profile has been greatly improved. 

To request high S/N of the stacked profile, we require the member galaxies in each 
c-R^Q-M volume to be at least 30. With this constraint, we get 193 composite galaxies, 
containing ~ 80 percent of the initial 25,517 individual galaxies. Those missed galax- 
ies are located in the outskirts of the c-R^q-M space with very low number density. 



2.3 The color gradient of disk component 

As shown in Fig.[2l the color of this composite galaxy becomes bluer with increasing 
radius. The color profile is steeper in the central part and becomes shallower outwards. 
The outer part of color profile (R > R^q) shows a nearly linear shape. 



Following the study of iTaylor et al.l (|2005[) . we separate the color profile of each 
composite galaxy into two parts at R^o and consider the region R > R^q as the disk 
component. We measure the color gradient of disk component within the region from 
i?50 to 3R50. The outer boundary is based on the considerations that a star formation 
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Fig. 2 The g — r (upper panel) and r — z (lower panel) color profiles of an 
example composite galaxy which is composed of 166 member galaxies with 
1.9 < c < 2.0, 7.0 < R50 < 8.0 kpc and -21.5 < M < -21.0 mag. 
Each small dot denotes the color at the radius of photometric annulus for a 
member galaxy. The dashed line is the best linear fitting of all the small dots 
between R^q ~ 7.5 kpc and SR^q (denoted by two vertical dotted lines). The 
large circles show the average color at given radius bins, with open and filled 
circles separated at R50. The error-bars show the scatters. 



threshold occurs outside the radius of about 5-tim es of the scale-length R^ (~ S-Rso) 
(Ivan der Kruit & Searld . Il98ll : IPohlen et al.L |2000|) and the large uncertainties of the 
photometric measurements exist in the very outer regionj^. 



The typical errors of surface brightness at 3i?5o are 0.105 for gr-band, 0.095 for r-band and 0.273 
for z-band, respectively. 
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Since the color profile of disk component (R^q ~ 3i?5o) can be well parameterized 
by a straight line, we use the linear least squares fitting to model the gradients of 
galactic disks. In specific, we make a linear fitting of all the data points in the range 
-R50 ~ 3-R50 and estimate the fitting error with 20 bootstrap samples. The best linear 
fitting for the example galaxy is shown as the dashed line in Fig.[2l 

The disk color gradient is then defined as the slope of the fitting line, 



G 



kpc 



AC/ 
Ai?(kpc) 



(1) 



where CJ is the color index. For our composite galaxy, this Gkpc can be easily 
converted to the scaled color gradient, e.g. G50 (the color variation per -R50), by 
G50 = Gkpc X R50, where R50 is the average size of member galaxies in each sub- 
set. 

In this paper, we will only present the results of colors g — r and r — z. We exclude 
the colors related to u band because of its low image sensitivity and high sensitivities to 
recent star formation and dust attenuation. We do not use z-band because of the 'Red 
halo' effect of the Point Spr ead Funct i on (PS F), which was introduced by 



Michard 



(|2002|) and was first found by lWu et al.1 (|2005|) in SDSS z-band. E xcept for the i-b and, 
the PSF profiles in g, r and z bands match well with each other (|Wu et al.L l2005h . so 
the effects of the PSF mismatch on the g — r and r — z color profiles are negligible. 



3 RESULTS 



3.1 Color gradient 



We measure the color gradients of disks (R50 ~ 3.0i?5o) of all 193 composite galaxies 
and show their distributions in Fig. |3l Panels (a) and (b) show the histograms of g — r 
gradient Ggr and r — z gradient Grz, whereas panels (c) and (d) show the histograms 
of their errors, respectively. 



Consistent with many early studies (e.g.. Ide Jong 



Moth & ElstonI l2002l : iTaylor et all l2005l : iLi et al 



2005 



1996b 



Bell & de Jond l2000l: 



Bakos et al 



2008D, most of 

our composite galaxies show negative color gradients. We notice that, forg — r color, 
about 25 percent ones have positive gradients. As we will show below, these composite 
galaxies with positive Ggr are mainly those with small -R50 and large uncertainties. 

As can be seen in Fig. [3l both Ggr and Grz show roughly Gaussian distributions 
with means and scatters G 



gr 



-0.006 mag kpc~\ 



-0.018 mag kpc ^ and 



a 



gr 



0.009, arz = 0.011. The typical errors of fitted gradients are about 0.003 for 
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Fig. 3 The distributions of color gradients Ggr [panel (a)] and Grz [panel (b)] 
of composite galaxies. The bottom two panels show the distributions of errors 
of Ggr[ S Ggr, panel (c)] and Grzi S Grz, panel (d)] respectively. 



Ggr and 0.005 for Grz, which are much smaller than the scatters agr and arz- If we 
de-convolve these typical errors from the distributions of Ggr and Grz, we get intrin- 
sic scatters (Tgr,i = 0.008 and arz,i = 0.010. These non-negligible intrinsic scatters 
among different composite galaxies imply that the color gradient is intrinsically corre- 
lated with the physical properties of galaxies. We discuss possible correlations in detail 
below. 

3.2 Correlation with galaxy properties 

In this section, we study the correlations between the color gradients and structural 
parameters of galaxies and try to find out which one is the main factor regulating the 
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color gradients. Technically, we use the Pearson's correlation coefficient to quantify 
each possible correlatior0. In the case of our sample containing 193 data points, the 
confidence level of the correlation will be over 99 percent if the absolute value of a 
correlation coefficient is larger than 0.18. 

As described in section 12.21 our composite galaxies are constructed in c-R^q- 
M space, so we explore the possible correlations between color gradients and these 
three bin parameters first. Fig. |4] shows the measured color gradients (Ggr and Grz) 
against the r-band concentration c, physical size R^o and absolute magnitude M in 
the top, middle and bottom panels, respectively. Each circle in the figure represents 
one composite galaxy, with sizes and colors of the circles indicating R^q and M of 
the galaxies. Larger circles represent larger galaxies and redder colors indicate more 
luminous galaxies. The errors of the color gradients are not shown for clearness (but 
shown in Fig. [5]). The Pearson's correlation coefficients between color gradients and 
these three parameters (c, R^q and M) are labeled in the corner of each panel and listed 
in Table E 

Table 1 The Pearson's correlation coefficients between color gradients (Ggr 
and Grz) and four structural parameters (r-band): concentration parameter c, 
physical half-light radius i?5o, absolute magnitude M and disk surface bright- 
ness i^d- 





C 


R50 


M 




Ggr 


-0.08 


-0.24 (-0.08*) 


-0.09 (-0.33*) 


-0.71 


Grz 


0.01 


0.03 (0.05*) 


-0.39 (-0.46*) 


-0.75 



Notes: *- The correlation coefficients of composite galaxies with i?5o > 3kpc. 



3.2. 1 With concentration parameter c 



As we can see from the top two panels of Fig. 4, the color gradients are almost inde- 
pendent of the concentration parameter. The correlation coefficients rcgr,c = —0.08 
and rGr^,c = 0.01 confirm that the color gradients are not correlated with c. 

It is well known that the concentration parameter is rel ated to the morphologi 



cal type, with earlier type galaxies h a ving l arger values of c dShimasaku et al. . 



2001 



Strauss et al. . 



2002 



Nakamura et al. 



20031). For spiral galaxies, a larger value of c 



^ During the calculation of this coefficient, the error of the color gradient of each composite galaxy 
has been used as a weight. 
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Fig. 4 The color gradients {Ggr and Grz) as functions of concentration param- 
eter c, physical half-light radius R^q and absolute magnitude M. The corre- 
lation coefficients r are labeled in the lower-right comer of each panel. Each 
point represents the color gradient of one composite galaxy. Bigger circles 
denote galaxies with larger sizes and redder colors represent higher luminosi- 
ties. 



also indicates an earlier type spiral, i.e. larger contribution from bulge component. 
Therefore, our results imply that the color gradient of the disk component is not re- 
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lated to the bulge. Moreover, this independence of the bulge component also confirms 
that our separation of disk component at R^q is adequate to avoid the bulge contami- 
nation. 

3.2.2 With physical size R^o 

For physical size R^q, as can be noticed from the correlation coefficients, the gradient 
Ggr shows a weak correlation {rGg^,R^o — ~0.24), while Grz is independent of -R50 
('"Grz.Rso = 0.03). If we go to the plot of Ggr vs R50 [panel (c) of Fig. IH, we can 
see that the weak correlation between Ggr and R^o is mainly caused by the galaxies in 
the smallest size bin, where more than half of the composite galaxies show positive 
gradients. Actually, if we exclude these galaxies from the correlation analysis, we get 
a coefficient r = —0.008, a similar independent behavior as that of the r — z gra- 
dient. Meanwhile, the correlation coefficient rGr~,B.r,o changes from 0.03 to 0.05. We 
also mention that the measurement of the color gradients of these smallest galaxies 
is quite uncertain (see the error-bars of these smallest galaxies in Fig. [S]) because of 
the limited number of member galaxies and pro f Mean data points. However, the fact 
that the smallest galaxies show positive Ggr but negative G^,.^ may relate to the migra 



tions of star-burst regio ns in blue compact galaxies (e.g. JSchulte-Ladbeck et al.lll999 



de Paz & Madore 



2005|) and worths further studying. 
Another interesting point we can see from the figure is that, the dispersion of 
the color gradients (for both Ggr and Grz) tends to be smaller for bigger galax- 
ies. For the biggest galaxy {R^q ~ 9 kpc), the gradients roughly keep as constants 
Ggr.gkpc ~ —0.008 and Gjz^gkpc ~ —0.016 (where the dispersion is compa r able t o 



measurement error). Similar result has also been found by lMunoz-Mareos et al.l(|2007|) . 
who investigated 161 nearby face-on spiral galaxies and found smaller FUV — K gra- 
dient dispersion for bigger galaxies. 

3.2.3 With absolute magnitude M 

The correlation coefficients between the color gradients and absolute magnitude are 
^Ggr,M = —0.09 and rGr^,M = —0.39. Again, when we exclude the composite galax- 
ies in smallest size bin as above, the coefficient r^^^, a/ increases to —0.33 and rG^,,M 
becomes —0.46. The weak anti-correlation between the color gradient and absolute 
magnitude indicates that brighter galaxies have shallo wer gradients. This result is ap- 



parently conflict with the result of iKim & AnnI (|l990h . where they found more lumi 



nous galaxies tend to have steeper color gradients. However, this apparent contradic- 
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tory is caused by different definitions of the color gradients. Kim & AnnI (|1990|) defined 
the color gradient as color variation per scaled radius -02^ while we adopt the defi- 
nition of color variation per kpc. As we have mentioned, the color gradient Gkpc can 
be converted to the scaled color gradient by m ultiplying a factor of physical radius. 
Since more luminous galaxies have larger sizes (IShen et al.Ll2003|) . if we convert Gkpc 
to Gso, we also get a positive correlation between G50 and absolute magnitude with co- 
efficients 0.17 and 0.18 for g — r and r — z gradient respectively, i.e. the more luminous 
galaxies have steeper scaled color gradients. 

Till now, we have explored the correlations between the color gradients and three 
bin parameters separately and found a tendency of lower scatter of color gradients 
for larger galaxies and a weak correlation between the color gradient and absolute 
magnitude. If we combine these bin parameters together, we may find more interesting 
results. In the middle two panels of Fig. HI if we focus on symbol colors, we will find 
a systematical change that less luminous galaxies (with bluer symbol colors) tend to 
have steeper color gradients in a given radius bin. Similar systematical change can 
also be found in bottom panels of Fig. HI where larger galaxies (with bigger symbol 
sizes) have steeper color gradients at given absolute magnitude. These two phenomena 
indicate a strong correlation between the color gradients and surface brightness. We 
discuss this correlation in detail below. 



3.2.4 With disk surface brightness 

In this section, we turn to investigate the correlation between the color gradients and 
disk surface brightness fid- Here, is defined as average surface brightness from R^q 
to Rqq and has been corrected by K-correction and cosmological dimming effect. This 
definition of the disk surface brightness avoids possible contaminations from bulge 
component and is consistent with our measurement of the disk color gradient. We 
calculate fid for each member galaxy and then average them to get the disk surface 
brightness of the composite galaxy. 

We plot the color gradients against disk surface brightness in Fig. \5\ Both Ggr and 
Grz show a remarkable anti-correlation with fid, with correlation coefficients rGgr,ii.a = 
—0.71 and tq^-^^i^ = —0.75. We use the least squares method to fit linear relations 
between the color gradients and fid weighted with the errors of the color gradients. The 



Kim & And (Il990h defined the -D25 as isophotal diameter for /i = 25magarcsec ^ isophote in 



B-band. 
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Fig. 5 The relationship between color gradients {Ggr and Grz) and disk sur- 
face brightness /i^. Each point represents the color gradient of one composite 
galaxy and the error-bar shows the measurement error estimated from boot- 
strap samples. The solid line in each panel shows the best linear fitting of all 
the data points and the dotted lines present the region of ±1(t rms scatter. 



best fitting results are 

Gar = -0.011Ud + 0.233, 

(2) 

Grz = -0.015/irf + 0.324, 

which are shown as the solid lines in Fig. \5\ The rms{a) of these fitted relations are 
0.005 for Ggr and 0.006 for Grz- The range of the ±lcr region is shown as two dotted 
lines in each panel. As we can see, most of the galaxies outside the ±la region are 
small galaxies (small symbol size) with large measurement errors. Moreover, we have 
tested that the residuals of these best fitting relations are not correlated with either 
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the luminosities or the sizes of galaxies (see the colors and sizes of symbols). That 
means, the weak correlation between the color gradients and absolute magnitude we 



have found in Section 13.2.31 is actually caused by this G — fid relation and known 



L — fid relation (e.g., iShen et al.l 120031 ). This also indicate that the tight correlation 
between color gradients and disk surface brightness is an intrinsic relation of spiral 
galaxies, which gives a strong observational constraint on the models of formation and 
evolution of galactic disk. 

4 SUMMARY AND DISCUSSION 

In this study, a sample of 25,517 face-on spiral galaxies drown from SDSS-DR4 is used 
to investigate the color gradients of galactic disks. We construct composite galaxies to 
improve the S/N of color profiles and measure their color gradients, and then investi- 
gate the correlations between the color gradients and structural parameters, including 
concentration, physical size, luminosity and disk surface brightness. The main results 
are listed as follows. 

1 . The color profile of disk component, from R^q to S-Rso, can be well represented by 
a straight line. Almost all of composite galaxies show negative r — z gradient, while 
about 75 percent of them show negative g — r gradients. The mean color gradients 
are Ggr = —0.006 mag kpc^^ and Grz = —0.018 mag kpc~^, respectively. 

2. The color gradients of galactic disks are independent of morphological types (con- 
centration c). 

3. There is also no significant correlation between disk color gradients and the size of 
galaxies (-R50), but bigger galaxies tend to have smaller scatters and uncertainties 
of the color gradients. 

4. The color gradients are dependent on the absolute magnitude M, with fainter galac- 
tic disks having steeper color gradients. But this correlation may caused hy G — fXd 
relation. 

5. The most remarkable correlation we find is between the color gradients and disk 
surface brightness fid, which can be fitted by straight lines: Ggr = —O.Ollfid + 
0.233 and Grz = -O.OlS/i^ + 0.324. 

There are two potential causes of color gradients in galactic disks. One is the sys- 
tematical change of dust attenuation along the disk and the other is the gradient of 
stellar population. 
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For dust, a reasonable and simplified assumption is that the du st distribution along 
the disk is proportional to the stellar mass density profile (e.g., IPhillipps & Disney 

Thus, the inner disk contains more dust than outer region, and 



1988 



Regan et al 



naturally generate a negative color gradient (redder inward) as we observed. However, 
in this scenario, high surface brightness (surface mass density) galaxies have more dust 
and will generate steeper color gradients than low surface brightness galaxies, which 
is opposite to our observational results. Therefore, dust attenuation could not be the 
only role in the origin of color gradients of galactic disks. Many previous works also 
have suggested that the dust only plays a minor role in regul a ting the color gradients 
of spiral galaxies (e.g., Ide Jonglll996d : iBell & de Jongll2000l : 



MacArthur et al. 



2004: 



Taylor etal.ll2005 ). 



Regarding the stellar population, it is related to chemical evolution history of the 



spiral galaxies arid has only been studied in details for 



MW: 



Chang et al 



199 



Chen& 



ration; M33: iMagrini et al. 



2007 



Hou 



2003 



Fu et al 



2009 



ocal group galaxies (e.g., 
; M31: Yin et al. in prepa- 



Barker & Sarajedinil l2008h . It is broadly accepted 



that the star formation rate (SFR) is propo rtion to the surface mass density (e.g.. 



Bell & de 



Jong 



2000 



Kennicutt 



Kauffmann et al, 



19981: 



2003b. and rn av also inverse proportion to ra- 



Boissier & Prantzos 



19991), which means SFE (star forma- 



dius (e.g., 

tion efficiency, defined as SFR/Eg^^, where T,gas is gas surface density) is higher in 
the inner region than outer part. Thus, in the inside-out star-forming scenario, the inner 
regions of galactic disks evolved faster than outer regions due to higher SFE, leading 
to older and mor e metal-rich stellar populat ions and generating negative color gradi- 



ents as observed (Prantzos & Boissier 



2000L Yin et al. in preparation). Moreover, the 



Kennicutt st ar formation low is nonlinearly proportional to the surface brightness as 
SFR oc fi^-^ (|Kennicuttlll998h . So, the higher surface bri ghtness galaxies have higher 
SFE, i.e. the chemical evolution will be more adequate. iBoissier & PrantzosI (|2000|) 
have modeled the color profiles of galactic disks using a nonlinear star formation law 
and found that more massive and smaller galaxies (i.e., higher surface density ) favor 
a more rapid early star formation than the lower surface density ones. Their gas frac- 
tions are lower and chemical evolutions are more adequate nowadays even in the outer 
region, which means the abundance and stellar population vary with radius slowly. 
Therefore the color differences between the inner and outer regions of high surface 
brightness galaxies will be smaller than that of low surface brightness ones. This sim- 
ple argument concludes a consistent trend with our observational results that the higher 
surface brightness galaxies have shallower color gradients. However, besides this static 
chemical evolution, the growth and evolution of galactic disks are related to a lot of 
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complicated processes , e.g. gas cooling, outflow, radial mixing, minor merging (e.g.. 



Kaiser & Binneyll2003l : iBlock et al.ll2006l) . To fully understand and quantify the origin 
and evolution of color gradients observed in spiral galaxies, more detailed models are 
required. 
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